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ABSTRACT 

The  study  represents  a  continuation  of  the  development  of 
a  model  output  statistics  scheme  to  specify  marine  fog  over 
the  open  ocean  and  in  coastal  waters.   Thirty-seven  direct 
and  derived  Fleet  Numerical  Weather  Central  model  output  para- 
meters, monthly  climatological  fog  frequencies,  combinations 
of  the  aforementioned  parameters  (i.e.,  interactive  parameters) 
and  a  persistence  parameter  are  used  as  predictors  in  a  step- 
wise multiple  linear  regression  approach  to  estimate  a  predic- 
tand  defined  as  marine  fog  probability.   The  predictand  is 
categorized  in  two  ways,  in  one  case  (FOGCAT  I )  as  smoothed 
probabilities  from  0  to  100%  as  a  function  of  present  weather, 
past  weather,  visibility  and  low  cloud  type;  and,  in  another 
case  (FOGCAT  IT)  as  a  limited  number  of  discrete  probabilities 
(to  include  0  and  100%)  derived  from  present  weather,  past 
weather  and  visibility  only.   This  study  derives  diagnostic 
regression  equations  only  using  as  a  dependent  data  sample 
over  24,000  North  Pacific  Ocean  (30-60N)  surface  synoptic  ship 
observations  at  0000  GMT  for  June  through  August  1976  and  1977. 
The  predictor  parameters  contributing  most  significantly  to 
the  variance  are  sensible  and  evaporative  heat  fluxes,  monthly 
climatological  fog  frequencies,  and  meridional  wind  speed. 
Threat,  Heidke  skill,  and  Panof sky-Brier  probability  scoring 
methods  are  applied  to  a  selection  of  the  derived  equations. 
Predictand  variance  explained  reaches  .670,  threat/skill 
scores  reach  0.42  and  probability  scores  are  as  low  as  0.28 
using  the  FOGCAT  I  categorization  scheme  for  the  predictand. 
Equations  for  June  and  July  appear  more  stable  than  those  for 
August . 
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I.   INTRODUCTION 

Marine  fog  presents  a  threat  to  safe  nautical  activities 
as  well  as  to  low-level  aviation  over  the  oceans.   Histori- 
cal monetary  and  human  losses  associated  with  United  States 
Navy  operations,  attributable  solely  or  mostly  to  fog,  have 
been  documented  by  Wheeler  and  Leipper  (1974).   These  types 
of  losses  should  be  significantly  reduced  as  the  methods  for 
analysis  and  forecasting  of  marine  fog  become  more  accurate. 
The  environmental  science  community  would  additionally  bene- 
fit from  the  increased  accuracy  of  fog  analyses/ forecasts  by 
a  commensurate  improvement  in  specifying  fog-associated  para- 
meters in  various  analysis  models,  especially  boundary  layer 
models . 

Research  in  the  past  several  years  conducted  by  the  De- 
partments of  Meteorology  and  Oceanography,  Naval  Postgraduate 
School  CNPS),  Monterey,  California,  has  been  directed  at  more 
adequate  diagnoses  of  marine  fog  and  a  better  representation 
of  its  climatology,  as  based  on  ocean-station-vessel   and 
ship-of -opportunity  observations  (Renard,  Englebretson  and 
Daughenbaugh,  1975;  Renard,  1976).   Addit ionally,  more  recent 
research  has  attempted  to  use  geostationary  weather  satellite 
data  to  diagnose  marine  fog  areas  (Ihli  and  Renard,  1977; 
McNab,  1979).   The  totality  and  continuity  of  such  satellite 
observations  makes  this  approach  potentially  very  useful. 


11 


However,  much  further  testing  is  needed  before  operational 
application  is  feasible. 

In  a  more  statistical  vein,  an  earlier  study  by  Nelson 
(1972)  developed  regression  equations  using  the  "perfect 
prognostic  approach"  to  forecast  visibility  (and  hence  fog) 
at  ocean  station  vessels  in  the  North  Atlantic  Ocean.   This 
approach  is  not  unlike  the  Model  Output  Statistics  (MOS) 
methods  used  currently  by  the  National  Weather  Service  (NWS) 
to  forecast  weather  parameters  of  interest  over  the  conti- 
nental U.S.,  Hawaii  and  Alaska  (Glahn  and  Lowry ,  1972),  ex- 
cept the  "imperfect  prognostic  approach"  is  the  present  mode. 
Nelson  found  quite  high  correlations  of  visibility  with  wind 
speed,  relative  moisture  content  and  various  evaporative 
parameters,  using  observed  data.   However,  when  he  attempted 
to  forecast  fog  using  output  parameters  from  the  United 
States  Navy  Fleet  Numerical  Weather  Central's  (FNWC)  numeri- 
cal prediction  models,  his  regression  equations  showed  little 
skill.   Earlier  work,  conducted  by  Schramm  (1966),  found  high 
correlations  between  observed  values  of  humidity,  air  tempera- 
ture, wind  speed  and  visibility  at  sea;  however,  the  regres- 
sion equation  developed  from  these  parameters  was  of  little  value 
at  times  of  low  visibility. 

Additional  research  into  the  feasibility  of  developing 
regression  equations  to  forecast  marine  fog,  based  on  FNWC ' s 
model  output  parameters,  has  been  conducted  by  Van  Orman  and 
Renard  (1977)  and  Quinn  (1978).   Both  studies  developed 
schemes  which  statistically  processed  up  to  37  direct  or 
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derived  FNWC  model  output  analysis-time   parameters  in  con- 
junction with  surface  synoptic  ship  observations  to 
generate  multiple  linear  regression  equations  yielding  as  a 

predictand  the  probability  of  marine  fog  occurrence  at 
analysis  time  over  the  North  Pacific  Ocean.   The  Van  Orman 
/Renard  results,  based  solely  on  July  1976  data,  showed 
considerable  skill  over  FNWC ' s  advective  fog  model,  FTER, 
as  applied  to  the  analysis  of  fog  (U.S.  Naval  Weather  Service 
Command,  1975)  and  climatological  fog  frequencies  developed 
at  the  NPS  CWillms,  1975).   Quinn  expanded  the  data  base  to 
also  include  June  and  August  1976  and  developed  regression 
equations  both  for  the  unique  probabilistic  marine  fog  pre- 
dictand introduced  by  Van  Orman  and  Renard  as  well  as  a 
modification  of  that  predictand.   Quinn  further  introduced 
interactive  predictor  parameters  into  the  regression  scheme 
by  combining  the  climatological  fog  frequency  parameter  with 
the  two  most  significant  model  output  parameters  to  derive  a 
new  set  of  diagnostic  regression  equations.   A  clear  improve- 
ment over  climatology  and  FNWC ' s  FTER  was  demonstrated  using 
these  modified  regression  equations  on  dependent  data.   The 
equations  based  on  the  June  1976  data  also  showed  improvement 
over  FTER  and  climatology  when  used  on  an  independent  data 
set  from  July  1976.   These  prior  studies  give  clear  indication 
of  the  relative  accuracy  of  the  MOS  approach  in  diagnosing 
marine  fog,  with  a  successful  extension  to  forecast  modes  yet 
to  be  demonstrated. 


Some  parameters  are  taken  from  analysis  models,  other  are 
diagnostic  parameters  from  prognostic  models. 
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Currently,  only  FNWC  is  producing  fog  analysis/forecasts 
on  a  large  scale  over  the  oceanic  areas  through  their  sta- 
tistical probabilistic  product,  FTER.   A  climatological 
parameter  is  not  used  in  their  product  due  to  the  unavail- 
ability of  an  accurate  climatology  at  the  time  of  its 
development.   A  comprehensive  uniform  climatology  of  fog 
occurrence  over  all  of  the  Northern  Hemisphere  is  now  avail- 
able from  the  National  Climatic  Center  (Guttman,  1978)  and 
has  been  incorporated  into  this  study.   This  inclusion  re- 
presents a  refinement  of  the  fog  parameter  from  the  two 
prior  studies . 
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II.   OBJECTIVES  AND  APPROACH 

The  primary  objective  of  this  study  was  to  continue  the 
development  of  a  multiple  linear  regression  approach  begun 
by  Van  Orman  and  Renard  (1977)  and  pursued  by  Quinn  (1978) 
to  specify  effectively  the  distribution  and  likelihood  of 
marine  fog  at  analysis  time  over  the  middle  latitude  (30°- 

60  N)  region  of  the  North  Pacific  Ocean  during  the  summer 
season.   Part  of  the  primary  objective  was  to  evaluate  the 
equations'  skill  in  comparison  to  the  climatological  fog  fre- 
quencies derived  by  Guttman  (1978)  and  FNWC ' s  existing  opera- 
tional scheme  for  forecasting  advection  fog  probability, 

FTER.   A  second  objective  was  to  determine  whether  the  pre- 
dictand  categorization  scheme  FOGCAT  I  developed  by  Van  Orman 
and  Renard  or  the  FOGCAT  II  scheme  developed  by  Quinn  is 
superior  in  defining  the  occurrence  of  fog.   Additionally, 
an  investigation  of  the  usefulness  of  persistence  as  a 
diagnostic  prediction  parameter  was  initiated.   The  approach 
follows  that  of  Quinn 's  (1978)  experiments  with  the  addition 
of  June,  July  and  August  1977  data. 

Sets  of  regression  equations  were  developed  for  both  fog 
categorization  schemes  for  (.a)  each  of  the  six  months  June,  July 
and  August  1976  and  1977;  (b)  the  combined  data  for  both 
Junes,  both  Julys  and  both  Augusts,  (c)  the  combined  data 
for  June,  July  and  August  1976,  (d)  the  combined  data  for 
June,  July,  and  August  1977,  and  (e)  all  of  the  data  for  the 
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six-month  period.   Those  equations  which  were  based  on  data 
from  both  June's,  both  July's,  both  August's,  and  all  six 
months  were  scored  using  three  methods,  skill,  threat  and 
probability . 
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III.   DATA 

A .  AREA 

The  area  of  study  was  confined  to  the  North  Pacific  Ocean 
to  maintain  continuity  and  consistency  with  the  previous 
studies  of  Van  Orman  and  Renard  (1977)  and  Quinn  (1978), 
specifically  the  area  north  of  30°N  and  south  of  60°N  (Figure 
1).   A  27x14  grid  was  superimposed  on  this  region  with  grid 
points  coincident  with  their  counterparts  on  the  same  subarea 
of  the  standard  FNWC  63x63  grid  (Figure  2). 

B.  TIME  PERIOD 

Due  to  the  high  frequency  of  marine  fog  in  the  summer 
months  over  the  North  Pacific  Ocean  (Willms,  1975;  Guttman , 
1978),  the  additional  summer  months  (June,  July  and  August) 

of  1977  were  added  to  the  available  set  of  the  same  months 
in  1976.   At  0000  GMT  the  entire  North  Pacific  Ocean  is  in 
daylight  with  local  noon  occurring  at  the  international  date- 
line.  Only  data  from  this  time  were  used  since  the  accuracy 
of  visibility  reports  from  ships  transitting  the  area  should 
be  best  during  daylight  hours , 

C.  SYNOPTIC  WEATHER  REPORTS 

The  June-August  1976  and  1977  synoptic  weather  reports 
used  in  this  study  were  provided  by  the  Naval  Weather  Service 
Detachment,  physically  located  with  the  National  Oceanic  and 
Atmospheric  Administration's  National  Climatic  Center, 
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Asheville,  North.  Carolina.   These  data  had  been  edited  to 
eliminate  duplicate  reports.   The  total  number  of  observa- 
tions received  were  4277  for  June  1976,  4391  for  July  1976, 

4134  for  August  1976,  4232  for  June  1977,  4198  for  July 
1977  and  3156  for  August  1977.   All  observations  were  ob- 
tained from  transitting  ships  and  Ocean  Weather  Station 
Papa  except  data  obtained  from  11  land  stations  located  in 
the  study  area  (Appendix  A).   As  all  of  these  land  stations 
are  located  on  islands  or  immediately  adjacent  to  the  coast 
and  have  relatively  low  elevations,  the  effect  of  topography 
is  thought  to  be  minimum. 

D.  MODEL  OUTPUT  PARAMETERS 

FNWC  provided  the  22  diagnostic  model  output  parameters 
(MOP's)  for  the  time  period  and  area  of  interest  as  output 
from  several  of  their  numerical  analysis  and  prediction 
models:  the  Mass  Structure  Analysis  Model,  the  Primitive 
Equation  Prediction  Model,  the  Marine  Wind  Model,  and  the 
Spectral  Ocean  Wave  Model  (U..S.  Naval  Weather  Service  Command, 
1975).   Additionally,  a  set  of  15  parameters  were  derived 
from  the  primary  set  and  eight  interactive  parameters  were 
developed.   A  complete  listing  and  description  of  all  these 
parameters  is  located  in  Appendix  B .. 

E.  FOG  FREQUENCY  CLIMATOLOGY 

Previous  research  used  the  Naval  Postgraduate  School  "s 
North  Pacific  Ocean  marine  fog  climatology  (Willms,  1975). 
As  an  updated  hemispheric  fog  frequency  climatology  (Guttman, 
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1978)  is  now  available,  it  was  used  exclusively  throughout 
this  research.   Specific  values  for  fog  frequency  were  inter- 
polated to  an  89x45  grid  in  the  area  of  study,  which  is 
approximately  four  times  as  dense  as  the  FNWC  grid  over  the 
region.   This  finer  resolution  allowed  a  more  accurate  speci- 
fication of  the  fog  frequency  climatology.  The  fog  frequency 
climatologies  for  each  month  are  presented  in  Figures  3-5. 
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IV.   PROCEDURE 

A.   PREDICTAND  CATEGORIZATION  METHODS 

The  two  methods  of  predictand  (fog)  categorization  used 
throughout  this  study  are  based  on  four  elements  of  the 
synoptic  reports,  namely:  present  weather  (ww) ,  past  weather 

(W)  ,  visibility  (VIS),  and  low  cloud  type  (CL) .   Observation 
elements  other  than  present  weather  are  used  due  to  incon- 
sistencies in  the  reported  observations,  which  result  in 
part  from  current  rules  for  coding  observed  data  (U.S.  De- 
partments of  Commerce,  Defense  and  Transportation,  1969) 

CSee  Appendix  C).   For  example,  fog  is  not  coded  as  present 
weather  whenever  any  form  of  precipitation  is  occurring 
simultaneously.   Both  methods  of  fog  categorization  assign 
a  probability  of  fog  occurrence  to  each  synoptic  report, 
differing  mainly  by  the  number  and  combination  of  observa- 
tion elements  used  to  assign  the  probability  value. 

The  Fog  Categorization  I  CFOGCAT  I)  scheme  was  developed 
by  Van  Orman  and  Renard  (1977)  and  used  all  of  the  four 
elements  listed  above.   This  scheme  assigns  one  of  five 
major  fog  categories  and  one  of  three  subcategories  to  each 
synoptic  report  as  a  function  of  ww,  W,  and  CL.   Then  a  fog 
probability,  0  to  100%  is  assigned  to  each  observation  de- 
pending upon  the  major  category  and  subcategory  assigned. 
Appendix  D  outlines  FOGCAT  I. 
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The  Fog  Categorization  II  (FOGCAT  II)  scheme  developed 
by  Quinn  (1978)  uses  only  present  weather,  past  weather  and 
visibility  to  assign  one  of  seven  discrete  fog  probabilities 

(0%,  10%,  15%,  35%,  60%,  85%  and  100%)  to  each  synoptic 
observation.   This  scheme  is  more  direct  in  its  method  of 
categorization  and  assigns  fog  probabilities  based  on  a 
filtering  technique.   First  present  weather  is  evaluated. 
If  present  weather  definitely  indicates  either  the  presence 

or  absence  of  fog,  then  a  probability  of  100%  or  0%  is 
assigned,  respectively.   In  either  of  these  cases,  both  past 
weather  and  visibility  are  ignored.   For  all  other  present 
weather  codes,  past  weather  and  visibility  are  considered 
in  the  assignment  of  a  fog  probability.   Past  weather  is 
evaluated  first.   If  the  possibility  of  fog  is  eliminated  by 
the  past  weather  code,  a  score  of  0%  is  assigned  and 
visibility  is  not  considered.   If  the  past  weather  code  in- 
dicates the  possibility  of  fog,  then  the  visibility  is  taken 
into  account  in  assigning  one  of  the  choices  for  an  inter- 
mediate fog  probability  (i.e.  10,  15,  35,  60,  80).   The 
assignment  of  these  intermediate  values  is  somewhat  arbitrary. 
Details  of  FOGCAT  II  are  presented  in  Appendix  E. 

Neither  the  FOGCAT  I  nor  FOGCAT  II  schemes  represent  an 
ideal  categorization  scheme  for  fog  probability.   FOGCAT  I, 

on  the  one  hand,  could  assign  an  observation  a  fog  probability 
as  low  as  62.1%  even  though  the  present  weather  code  definitely 
indicates  the  presence  of  fog,  or  it  could  assign  a  proba- 
bility as  high  as  24.1%  when  fog  appears  very  improbable  from 
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the  synoptic  report.   FOGCAT  II  was  developed  to  achieve  a 
better  fit  of  probability  to  occurrence  of  event.   Except 
for  the  "100"  and  "0"  values,  the  FOGCAT  II  probabilities, 
like  those  of  FOGCAT  I,  are  somewhat  arbitrary.   In  any  case 
the  schemes  were  accepted  as  developed  by  the  previous 
authors,  recognizing  that  some  adjustments  may  improve  the 
physical  representation  of  the  categorization  along  with  its 
verification . 

B.  PREDICTOR  INTERPOLATION 

The  interpolation  method  used  to  determine  values  of  the 
predictors  at  the  observation  points  is  a  natural  bicubic 
spline  curvilinear  interpolation  scheme.   The  locally  de- 
veloped program,  SPLIN,  is  available  at  the  NPS  W.  R.  Church 
Computer  Center  and  was  used  throughout  the  study. 

C.  REGRESSION  SCHEME 

A  stepwise  multiple  linear  regression  program,  called 
BMDP2R  (University  of  California,  1975)  was  chosen  as  the 
means  for  deriving  the  marine  fog  diagnostic  equations.   The 
polynomial  regression  scheme  also  availabe  in  the  BMD  com- 
puter program  series  was  not  used  due  to  its  inability  to 
handle  data  sets  of  over  1000  cases.   As  each  month  has 
approximately  4000  synoptic  observations  and  combined  data 
sets  had  at  least  7000  cases,  the  usefulness  of  the  polynomial 
regression  scheme  applied  to  smaller  segments  of  the  data  was 
unlikely  to  give  accurate  results.   A  non-linear,  non-poly- 
nomial scheme  also  available  in  the  BMD  series  uses  excessive 
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computer  memory  for  a  45-variable  regression,  and,  therefore, 
was  likewise  not  used. 

BMDP2R  computes  a  sequence  of  multiple  linear  equations 
in  a  stepwise  manner.   At  each  step  one  variable  is  either 
added  or  removed  from  the  previous  step's  equation  as  dependent 
on  the  F-to-enter  and  F-to-remove  criteria.   In  the  forward 
selection  procedure  the  dependent  predictor  inserted  in  the 

equation  is  the  one  with  the  highest  coefficient  of  partial 

2 

determination  (AR  ).   At  each  step  the  regression  procedure 

reevaluates  the  variables  already  in  the  equation,  and  may 
find  that  a  variable  important  at  an  earlier  stage  may  be 
less  important  at  a  later  stage  due  to  the  intercorrelation 

of  variables.   In  this  procedure,  before  a  variable  is  added, 
the  variable  already  in  the  equation  with  the  partial  F 
value  is  dropped  if  this  latter  value  is  less  than  the  maxi- 
mum F-to-remove  value  (Wesolowsky,  1976).   Selection  proceeds 
until  no  candidate  variables  qualify  or  until  all  independent 
variables  have  been  used. 

The  BMDP2R  program  permits  the  specification  of  the 
tolerance  level.   If  the  tolerance  is  near  zero,  the  variable 
being  considered  is  close  to  being  a  linear  combination  of 
the  variables  already  in  the  equation.   Variables  with  low 

tolerances  could  cause  computational  difficulties  in  the 
method  of  calculation  (Wesolowsky,  1976).   Since  tolerance  is 
defined  as  one  minus  the  coefficient  of  determination,  the 
value  of  .01  was  used  as  the  minimum  acceptable  tolerance 
criterion  (University  of  California,  1975). 
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The  approach  taken  in  this  study  was  to  treat  the  y-inter- 
cept  (Y-INTCP)  as  a  variable  and  to  use  the  BMDP2R  program's 
default  values  of  4.0  and  3.9  as  the  values  for  the  minimum 
F-to-enter  and  maximum  F-to-remove  criteria. 

D.   VERIFICATION  SCORING 

Two  types  of  verification  scores  are  used  in  the  study  to 
test  the  skill  of  different  regression  equations  in  describing 
the  distribution  of  fog  probabilities.   The  first  type, 
exemplified  by  the  Heidke  Skill  Score  (HSS)  and  the  Threat 

Score  (TS),  is  used  to  test  the  effectiveness  in  specifying 
discrete  occurrences  of  fog.   The  second  type  (Panof sky-Brier 
Probability  Score)  tests  the  accuracy  of  probabilities  in 
estimating  the  likelihood  of  fog  occurrence  (PS)  (Panofsky 
and  Brier,  1958) . 

The  formulae  used  for  the  three  scores  (HSS,  TS ,  and  PS) 
are  given  in  Appendix  F.   Currently  the  Threat  Score  is 
popularly  used  by  the  Techniques  Development  Laboratory, 
National  Weather  Service  (Bermowitz  and  Best,  1978;  Miller  and 
Best,  1978),  while  the  Heidke  Skill  Score  has  been  traditional- 
ly used  by  meteorologists.   The  Panof sky-Brier  Probability 
Score  is  suitable  for  evaluating  the  effectiveness  of  a 
predictor  given  in  terms  of  a  probability  (0-100%) . 
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V.   RESULTS 

A.   REGRESSION  EQUATIONS 

The  sets  of  equations  generated  for  the  individual  months 
of  June,  July  and  August  1976  and  1977  as  well  as  combina- 
tions thereof,  using  categorization  schemes  FOGCAT  I  and 
FOGCAT  II,  are  listed  in  Tables  I -XII.   They  are  shown  in 
stepwise  order  and  include  Y-intercept,  regression  coeffi- 
cients and  the  amount  of  variance  explained  with  the  inclu- 

2        2 
sion  or  removal  of  each  variable  (.R   and  AR  )  .   Only  six 

steps  are  given  in  each  table  for  purposes  of  brevity.   In 

2 
all  cases  no  two  successive  steps  added  a  cumulative  AR   of 

greater  than  .005  after  the  sixth  step  until  a  variable  was 
removed.   Recognizing  that  variables  which  explain  less  than 
1%  additional  variance  are  normally  not  included  in  opera- 
tional regression  equations,  six  steps  were  presented  so 
that  all  variables  having  a  possible  effect  on  further  devel- 
opment could  be  identified. 

Since  the  climatological  synoptic  regime  varies  from 
month-to-month  and  from  year-to-year  the  regression  equa- 
tions vary  accordingly.   Several  similarities  are  evident, 
however.   First,  in  all  cases  for  both  FOGCAT  I  and  FOGCAT  II 
the  first  model  output  parameter  to  be  entered  is  one  of  the 

heat  flux  terms,  either  EHF  or  SEHF .   Second,  either  CLIMO 
or  an  interactive  parameter  containing  CLIMO  enters  all  the 
equations  by  step  3  with  the  exception  of  the  equations 
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derived  from  August  1977  data  and  the  equation  for  all 
months  for  FOGCAT  I.   Finally,  the  meridional  component  of 
the  wind,  v,  or  an  interactive  parameter  containing  v  enters 
all  of  the  equations  no  later  than  step  4  with  the  exception 
of  the  equations  derived  using  FOGCAT  I  on  the  combined 
July  1976/77  data  or  FOGCAT  I  and  II,  on  the  August  1976 
and  FOGCAT  I  on  the  August  1977  data.   Additional  similari- 
ties are  evident  at  later  steps  but  are  not  considered 
significant  to  this  discussion. 

In  all  cases,  the  first  three  non-constant  variables 
entered  into  the  regression  equations  accounts  for  between 

71.9%  and  92.5%  of  the  summation  of  coefficients  of  partial 

2 

determination  (AR  )  determined  by  the  entire  set  of  model 

output  predictors . 

As  the  magnitudes  of  the  coefficients  of  determination 

2 

(R    1  vary  greatly  between  those  equations  derived  using  the 

two  different  categorization  schemes,  a  direct  comparison 
of  the  effectiveness  of  the  equations  based  on  FOGCAT  I 
versus  those  based  on  FOGCAT  II  is  difficult.   Although  Table 
XIII  indicates  that  the  EAR2  for  FOGCAT  II  is  generally 
greater  than  this  quantity  for  FOGCAT  I,  comparisons  are  more 
properly  made  using  standard  verification  scoring  methods, 
as  in  Section  B  following. 

The  eight  interactive  parameters  used  in  the  regression 
analysis  were  formed  by  either  multiplying  the  parameters 
found  to  be  most  significant  by  Quinn  (1978)  by  climatology 
or  by  multiplying  various  of  these  most  significant  parameters 
by  other  significant  parameters.   The  inclusion  of  these 
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interactive  parameters  allows  climatology  to  enter  the  scheme 

at  an  earlier  step.   Originally,  climatology  was  too  highly 
correlated  with  the  heat  flux  terms  to  enter  the  equations 
alone,  and  thus  it  was  rejected  by  the  tolerance  criterion. 

Having  identified  the  most  probable  meteorological  indi- 
cators of  fog  occurrence  by  regression  analysis,  the  feasi- 
bility of  adding  a  persistence  parameter  was  investigated. 
Due  to  lack  of  time,  this  investigation  was  limited  to  the 
July  1976  and  1977  data  sets.   The  persistence  parameter 
(PERS)  was  based  only  on  those  observations  which  had  at 
least  one  synoptic  report  within  a  one  degree  latitude/longi- 
tude box  of  the  reported  position  on  the  previous  day.   This 
requirement  reduced  the  data  set  for  July  1976  from  4391  to 
1368  reports  and  for  July  1977  from  4198  to  1591  reports. 

Regression  analysis,  with  PERS  included,  was  then  per- 
formed on  the  July  data  sets,  both  individually  and  in  unison, 

using  both  categorization  schemes.   The  sets  of  equations  for 
these  data  sets  are  listed  in  Tables  XIV-XVI .   Persistence 
displayed  a  high  correlation  with  EHF,  SEHF  and  CLIMO  and 
was  rejected  for  inclusion  in  the  equations  by  the  tolerance 
criterion  on  this  basis  until  the  seventh  step  at  the  earliest. 
The  similarities  among  equations  which  were  noted  earlier  are 
decreased  when  PERS  is  included  with  only  the  heat  flux  term 
retaining  its  former  importance  in  all  cases.   This  instability 
introduced  into  the  equations  may  be  only  the  result  of  the 
reduced  data  set;  the  investigation  of  persistence  as  a  fog 

indicator  should  not  be  abandoned  until  it  is  tested  on  a 
larger  scale,  both  alone  and  in  interactive  combinations. 
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Fields  of  selected  model  output  parameters,  FTER  and  re- 
gression-computed fog  probabilities  for  the  arbitrarily 
selected  observation  time,  0000  GMT  3  August  1976,  are  shown 
in  Figures  6-12.   Considerable  frontal  activity  in  the  mid- 
latitude  central  and  western  North  Pacific  Ocean  is  indicated 
by  the  sea-level  pressure  analysis  (Figure  6).   The  two  most 
important  fog-related  parameters,  EHF  and  SEHF ,  are  shown  in 
Figures  7  and  8.   Maximum  positive  values  at  the  subtropical 
latitudes  give  way  to  a  broad  zonal  band  of  negative  and  low 
positive  values  centered  between  45-55N.   Another  important 
regression  parameter,  the  algebraic  value  of  meridional  wind 
speed,  v,  is  displayed  in  Figure  9.   This  field  closely  re- 
lates to  the  geostrophic/gradient  wind  implied  by  the  sea- 
level  pressure  analysis.   Figures  10  and  11  show  the  regres- 
sion-computed fog  probabilities  based  on  three  variables  for 

Fog  Categorization  Schemes  I  and  II.   The  configuration  of 
these  fields  resembles  that  of  SEHF,  the  FNWC  variable  ex- 
plaining most  of  the  variance  in  fog  probability.   Probabili- 
ties on  the  FOGCAT  I  version  (Figure  10)  appear  to  be  about 
15%  greater  than  those  of  FOGCAT  II  (Figure  11).   The  differ- 
ence is  particularly  noticeable  in  the  subtropics  where  fog 
likelihood  is  climatologically  small.   The  FTER  probability 
field  CFigure  12)  shows  a  smaller  range  of  values  than  the 
NPS  probabilities,  and  is  a  closer  fit  to  the  FOGCAT  II 
regression  probabilities. 

B.   VERIFICATION  SCORING 

For  purposes  of  verifying  fog  occurrence  in  an  observation, 
major  categories  S  and  F  were  used  in  the  FOGCAT  I  scheme, 

28 


and  an  assigned  probability  of  fog  occurrence  of  100%  was  used 
in  the  FOGCAT  IT  scheme  CAppendices  D  and  E).   An  attempt  was 
made  to  use  85%  fog  likelihood  as  also  verifying  fog  occur- 
rence in  FOGCAT  II,  but  this  criterion  yielded  significant- 
ly poorer  results  than  100%  used  alone,  as  had  been  found 
by  Quinn  (1978).   This  may  indicate  that  some  fine  tuning  of 
the  assigned  fog  probabilities  in  FOGCAT  II  would  yield 
better  results. 

Optimum  threshold  probabilities  were  also  computed  for 
those  equations  based  on  the  combined  June  data,  combined 
July  data,  combined  August  data,  and  combined  data  for  June, 
July  and  August  1976  and  1977.   The  threshold  probability  is 

that  equation-specified  probability  which  best  separates  the 
fog/no  fog  occurrences,  hence  yielding  the  best  skill  and 
threat  scores  (Bermowitz  and  Best,  1978).   The  threshold 
Heidke  Skill  Scores  and  Threat  Scores  presented  in  Tables 
XV II -XX I  are  those  which  maximize  the  indicated  scores.   The 
Panof sky-Brier  Probability  Score  presented  in  these  tables 
directly  evaluates  the  accuracy  of  regression  probabilities 
in  estimating  the  likelihood  of  the  event. 

The  regression  equations  as  derived  by  the  BMDP2R  program 
were  not  used  in  their  entirety  in  the  verification  phase  of 
this  study  nor  were  all  the  equations  verified.   Only  the 
equations  in  the  form  containing  the  Y-intercept  and  the 
first  three  non-constant  variables  were  used.   It  is  possible 
to  "overfit"  the  regression  equation  to  the  dependent  data. 
Therefore  only  variables  which  bear  a  physical  relationship 


29 


to  the  predictand  should  be  included  (Panofsky  and  Brier, 
1958).   In  all  the  equations  verified,  the  first  three  vari- 
able predictors  included  had  an  apparently  strong  physical 
relationship  with  the  occurrence  of  fog,  while  several  equa- 
tions introduced  terras  with  only  a  weak  physical  relation- 
ship to  fog  likelihood  at  step  5. 

The  verification  scores  for  the  five  data  sets  for  which 
regression  equations  were  evaluated  and  compared  to  the  veri- 
fication scores  for  FTER  and  climatology  are  presented  in 
Tables  XVII-XXI.   Note  that  the  scores  for  the  regression 
equations  as  well  as  for  FTER  and  climatology  vary  from  the 
FOGCAT  I  to  the  FOGCAT  II  schemes.   Since  verification  in 
each  scheme  is  against  the  actual  occurrence  of  fog  (Category 
S  or  F  in  FOGCAT  I  or  100%  in  FOGCAT  II),  a  relative  compari- 
son of  the  effectiveness  in  categorizing  fog  occurrence 
between  FOGCAT  I  and  II  for  purposes  of  regression  analysis 
is  possible.   This  is  not  to  imply  that  observations  verify 
as  fog/no  fog  equally  with  FOGCAT  I  and  FOGCAT  II;  there  are 
slight  differences  not  considered  important  in  this  analysis. 

In  four  of  the  five  cases  verified,  the  skill  and  threat 
scores  for  the  regression  equations  were  higher  for  FOGCAT  I 
than  for  FOGCAT  II.   In  all  verification  scores,  the  FOGCAT 
II  regression  equation  outperforms  both  climatology  and  FTER, 
except  for  the  regression  equation  based  on  the  two  August 
data  sets.   The  FOGCAT  I  regression  equation  outperforms 
climatology  and  FTER  in  skill  and  threat  scores  for  each  data 
set  except  the  full  June,  July  and  August  1976  and  1977  set 
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where  climatology  is  the  top  performer.   In  probability  scor- 
ing, FTER  outperforms  the  FOGCAT  I  regression  equation  in 
four  of  the  five  cases  and  climatology  outperforms  it  in  the 
other  cases  in  the  P-score. 

The  variations  of  skill  and  threat  scores  for  climatology 
from  FOGCAT  I  to  FOGCAT  II  may  be  taken  as  a  measure  of  the 
differences  in  the  verifying  scheme  (i.e.  strong  fogger, 
category  S,  and  f oggers ,  category  F,  for  FOGCAT  I  and  observa- 
tions assigned  100%  probability  in  FOGCAT  II).   Thus,  a  mea- 
sure of  the  relative  comparison  of  FOGCAT  I  and  II  regression 

equations  might  be  obtained  by  noting  the  changes  in  the 
relative  improvement  over  climatology.   For  example  with 
reference  to  Table  XVII,  FOGCAT  I  skill/threat  scores  exceed 
climatologies  skill/threat  scores  by  .107/. 072,  while  for 
FOGCAT  II  the  numbers  are  only  . 082/ . 063, indicat ing  greater 
merit  for  FOGCAT  I.   Table  XXII  tabularizes  this  information. 
On  this  basis  relative  improvement  over  climatology  is  best 
using  FOGCAT  I  for  equations  in  Tables  XVII,  XVIII,  XIX 
(except  for  threat  score)  and  XXI,  while  FOGCAT  II  is  better 
for  equations  in  Table  XX. 

When  the  persistence  parameter  is  included  in  the  regres- 
sion scheme  for  the  two  Julys,  regression  skill  and  threat 
scores  generally  show  improvement  (Tables  XVIII-XIX).   Con- 
sidering the  reduced  data  set,  at  best,  it  is  to  be  regarded 
as  giving  encouragement  to  the  future  refinement  of  a  more 
comprehensive  persistence  parameter. 
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It  is  to  be  noted  again  that  all  regression  equations 
apply  to  analysis  time,  hence  are  diagnostic.   Further  all 
are  relative  to  dependent  data.   The  data  were  considered 
too  limited  to  expend  even  one  month  as  an  independent  set 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

From  the  evaluation  of  these  diagnostic  Model  Output  Sta- 
tistics regression  equations,  it  can  be  seen  that  the  MOS 
approach  to  prediction  of  marine  fog  holds  great  potential. 
The  developed  regression  equations  show  improvement  over  both 
climatology  and  the  FNWC  advective  fog  forecasting  product, 

FTER.   However,  it  remains  difficult  to  recommend,  without 
qualifications,  the  use  of  either  FOGCAT  I  or  FOGCAT  II 
schemes  for  the  predictand. 

The  comments  that  follow  relate  to  the  equation  for  each 
month  (for  two-year  period);  the  equation  for  all  six  months 
combined  are  not  considered  to  be  operationally  useful. 
Even  though  FOGCAT  II  equations  score  better  than  CLIMO  and 
FTER  (except  for  August)  using  all  the  scoring  methods, 
threat,  skill  and  probability,  FOGCAT  I  equations  show  higher 
threat  and  skill  scores  than  FOGCAT  II  and  excel  CLIMO  and 
FTER  in  two  of  the  scoring  methods,  threat  and  skill.   It 
appears  quite  significant  that  FOGCAT  I  shows  most  skill 
relative  to  climatology  and  FTER  as  well  and  therefore  this 
predictand  categorization  appears  to  be  the  best  to  use. 
Further  development  of  FOGCAT  I  and  II  or  another  approach 

may  well  yield  still  further  improvement  in  regression  results 

The  following  recommendations  are  offered  as  a  guide  to 
future  work  in  this  area: 
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1.  In  order  to  simplify  the  analysis  procedure  and  the 
amount  of  data  manipulation,  parameters  which  have  little  or 
no  effect  on  the  outcome  of  the  regression  equation  should 
be  deleted.   These  include:  SOLARAD ,  DDW ,  PPW ,  PDW,  SPW, 
SDW,  WCP,  THETAX,  THETAR,  STABX ,  ASTDX,  ASTDR,  CAPV,  ADTSEA , 
AASTDX,  AASTDR,  and  SSTA .   See  Appendix  B. 

2.  The  development  of  a  more  realistic  persistence  para- 
meter should  be  continued  and  included  in  the  regression 
scheme  both  alone  and  in  interactive  combinations.   As  a 
first  approximation  the  fog  probability  assigned  by  the 
regression  equation  (without  inclusion  of  a  persistence  para- 
meter) to  each  grid  point  of  the  previous  day  could  be  inter- 
polated to  the  actual  observation  positions  of  the  present 
day.   This  would  give  a  persistence  value  to  each  observation 
whether  or  not  a  report  was  available  within  a  one  degree 
latitude/longitude  box  on  the  previous  day.   These  first 
approximations  could  then  be  weighted  for  those  observations 
which  had  a  report  within  the  specified  box  on  the  previous 
day.   In  this  manner  the  size  of  the  data  set  would  not  be 
reduced  as  it  was  in  this  study, 

3.  Additional  data  (perhaps  one  more  year)  should  be 
added  to  further  stabilize  the  regression  analysis. 

4.  A  stratification  of  the  set  of  regression  equations 
by  latitude,  longitude,  or  by  meteorological  phenomena  (i.e. 
wind  direction,  positive/negative  thermal  advection,  or 
places  where  climatological  fog  frequencies  are  greater)  is 
likely  to  improve  accuracy. 


34 


5.  The  regression  approach  should  be  extended  to  regions 
other  than  the  North  Pacific  Ocean  and  time  periods  other 

than  summer. 

6.  As  FNWC  predictive  model  output  parameters  fields 
become  available  for  testing  they  should  be  incorporated  into 
the  regression  equations  for  forecast  intervals  to  at  least 
48  hours. 
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TABLE  XVIL   Skill,  threat  and  P-scores  June  1976  and  1977 
data. 


A.   FOGCAT  I 

Cate 

gorization 

Sc 

■heme . 

9 

R"  = 

.627, 

REGRESSION 
EQUATION 

FTER 

CLIMATOLOGY 

Skill  Score 

* 
.315 

.178 

.208 

Threshold 

40 

14 

27 

Threat  Score 

* 
.  338 

.258 

.266 

Threshold 

38 

21 

16 

P-Score 

.292 

.209 

* 
.178 

Fog  Probability  (%)  ■  36.725  -  2.095  (EHF) 

-  .030  (CLIMO)(SEHF)  +  .226  (CLIMO) 


B.   FOGCAT  II  Categorization  Scheme.   R  =  .367. 


REGRESSION 
EQUATION 

FTER 

CLIMATOLOGY 

Skill  Score 

.302 

.207 

.220 

Threshold 

28 

23 

27 

Threat  Score 

* 
.322 

.250 

.259 

Threshold 

26 

21 

22 

P-Score 

* 
.288 

.339 

.300 

Fog  Probability  (%)  -  31.154  -  2.956  (EHF) 

-  .032  (CLIMO)(SEHF)  +  .023  (CLIMO)(V) 


* 
Indicates  best  score  in  each  category 
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TABLE  XVIII.  Skill,  threat  and  P-scores  July  1976  and  1977 
data. 


A.   FOGCAT  I 

Cate 

gorization  Scheme, 

.  R2  = 

.670. 

REGRESSION 
EQUATION 

FTER 

CLIMATOLOGY 

Skill  Score 

* 
.384 

.336 

.303 

Threshold 

46 

35 

35 

Threat  Score 

* 
.404 

.355 

.360 

Threshold 

42 

23 

28 

P-Score 

.330 

* 
.216 

.242 

Fog  Probability  (%)  =  34.898  -  1.663  (EHF) 

-  .026  (CLIMO)(SEHF)  +  .292  (CLIMO) 


B.   FOGCAT  II  Categorization  Scheme.   R2  =  .441 


REGRESSION 
EQUATION 

FTER 

CLIMATOLOGY 

Skill  Score 

.378* 

.333 

.302 

Threshold 

38 

35 

42 

Threat  Score 

* 
.386 

.346 

.347 

Threshold 

34 

23 

33 

P-Score 

* 
.299 

.322 

.322 

Fog  Probability  (%)  =  37.131  -  4.037  (EHF) 

+  .029(jeLIMO)  (V)  +  .466  (EAIR) 


Indicates  best  score  in  each  category 
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TABLE  XIX.   Skill,  threat  and  P-scores  July  1976  and  1977 
data  with  persistence  parameter. 


A.   FOGCAT  I 

Cat' 

egorization 

Sch 

eme . 

R^  = 

.673, 

REGRESSION 
EQUATION 

FTER 

CLIMATOLOGY 

Skill  Score 

* 
.405 

.334 

.  319 

Threshold 

48 

28 

28 

Threat  Score 

* 
.423 

.371 

.391 

Threshold 

47 

25 

32 

P-Score 

.356 

* 
.237 

.281 

Fog  Probability  (%)  =  17.869  -  2.945  (EHF) 

+  .510  (CLIMO)  +  1.113  (EX) 


B.   FOGCAT  II  Categorization  Scheme.   R"  =  .458 


REGRESSION 
EQUATION 

FTER 

CLIMATOLOGY 

Skill  Score 

* 
.394 

.337 

.  320 

Threshold 

36 

33 

41 

Threat  Score 

* 
.414 

.365 

.  381 

Threshold 

36 

25 

33 

P-Score 

* 
.323 

.370 

.350 

Fog  Probability  (%)  =  28.167  -  3.956  (EHF) 

+  .642  (TX)  -  130.026  (STABR) 


Indicates  best  score  in  each  category. 
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TABLE  XX  . 


Skill,  threat  and  P-scores  August  1976  and 
1977  data. 


A.   FOGCAT  I 

Gate 

gorization  Sc 

:heme . 

IT  = 

.660 

REGRESSION 
EQUATION 

FTER 

CLIMATOLOGY 

Skill  Score 

* 
.403 

.374 

.146 

Threshold 

45 

28 

26 

Threat  Score 

* 
.378 

.  345 

.219 

Threshold 

40 

23 

14 

P-Score 

.279 

* 
.175 

.202 

Fog  Probability  (%)  -  33.781  -  1.859  (SEHF) 

+  1.059  (V)  -  .708  (CAPU) 


B.   FOGCAT  II 

Cate 

gorization  S 

cheme 

.   R   = 

.420 

. 

REGRESS  ION 
EQUATION 

FTER 

CLIMATOLOGY 

Skill  Score 

* 
.402 

.375 

.  136 

Threshold 

46 

28 

25 

Threat  Score 

* 
.  368 

.341 

.208 

Threshold 

40 

24 

16 

P-Score 

.297 

.280* 

.336 

Fog  Probability  (%)  =  26.352  -  2.105  (SEHF) 

+  .727  (V)  -  .067  (SEHF)(V) 


Indicates  best  score  in  each  category. 
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TABLE  XXI  .   Skill,  threat  and  P-scores  June,  July  and  August 
1976  and  1977  data. 


A.   FOGCAT  I 

Categ< 

orization 

Scheme . 

R   ■ 

.646, 

REGRESSION 
EQUATION 

FTER 

CLIMATOLOGY 

Skill  Score 

.214 

.196 

* 
.293 

Threshold 

27 

23 

44 

Threat  Score 

.267 

.251 

* 
.319 

Threshold 

16 

21 

41 

P-Score 

.304 

* 
.201 

.208 

Fog  Probability  (%)  =  31.966  -  1.559  (SEHF) 

+  .561  (V)  +  .215  (CLIMO) 


B.   FOGCAT  II  Categorization  Scheme.   R   =  .400. 


REGRESSION 
EQUATION 

FTER 

CLIMATOLOGY 

Skill  Score 

* 
.370 

.321 

.294 

Threshold 

32 

34 

40 

Threat  Score 

* 
.382 

.335 

.337 

Threshold 

31 

23 

32 

P-Score 

.283* 

.322 

.318 

Fog  Probability  (%)  =  30.406  -  2.811  (SEHF) 

+  .029  (CLIMO) (V)  +  1.530  (SHF) 


Indicates  best  score  in  each  category 
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TABLE  XXII: 


Comparison  of  Skill  and  Threat  Scores  for 
the  Regression  Equations  and  Climatology 


Period 

FOGCAT  I 

FOGCAT  II 

Skill  Score 
Difference 

(Regression 
Equation  - 
Climatology) 

Threat  Score 
Difference 

( Regression 
Equation  - 
Climatology) 

Skill  Score 
Difference 

(Regression 
Equation  - 
Climatology ) 

Threat  Score 
Difference 

(Regression 
Equation  - 
CI imatology ) 

June 
76  h   77 

.107 

.072 

.082 

.063 

July 
76  &   77 

.081 

.044 

.076 

.039 

July  76 
&  77  with 
persis- 
tence 

.086 

.032 

.074 

.03  3 

August 
76  &  77 

.257 

.159 

.266 

.160 
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APPENDIX  A 

LAND  STATIONS  USED  IN  STUDY 
(U.S.  Air  Force,  1972) 


Station 

Surface 

Number 

Latitude 

Longitude 

Elevation  (ft) 

32174 

45.2 

N 

147.9 

E 

38 

32186 

46.2 

N 

150.5 

E 

73 

32195 

46.9 

N 

151.9 

E 

26 

32207 

43.3 

N 

153.3 

E 

55 

32213 

50.9 

N 

156.7 

E 

42 

32215* 

50.7 

N 

156.2 

E 

unknown 

32217 

50.0 

N 

155.4 

E 

11 

32559 

53.1 

N 

160.0 

E 

88 

32594 

51.5 

N 

156.5 

E 

6 

32618 

54.9 

** 

N 

166.2 

** 
E 

19 

70454 

51.9 

N 

176.7 

W 

14 

Station  not  listed  in  reference.   Lat/long  values 
listed  are  as  given  with  each  synoptic  report 
furnished  by  ENWC. 


** 


Location  given  slightly  removed  from  actual  loca- 
tion of  55.2  N  latitude,  166.0  E  longitude  to 
accommodate  study  area. 


70 


APPENDIX  B 
OUTPUT  PARAMETER  DESCRIPTIONS 


Source  Model 


Symbol 


Name 
Description 


Units 


Analysis  Parameters  (FNWC's  Mass  Structure  Model) 

PS       Sea  Level  Pressure  (mb) 


Sea  Level  Pressure 

Analysis  of  observed  sea  level 

parameter. 


TAIR     Surface  Air  Temperature 

Analysis  of  observation-level 
air-temperature. 

EAIR     Surface  Vapor  Pressure 

Analysis  of  observation-level 
vapor  pressure  derived  from  the 
dew  point. 

T925     925  mb  Air  Temperature 

Analysis  of  925  mb  air  tempera- 
ture. 

TSEA     Sea  Surface  Temperature 

Once-daily  analysis  of  observed 
sea-surface  temperature. 


(°C) 


(mb) 


(°C) 


(°C) 


B.   P.E.  Parameters  (FNWC's  Primitive  Equation  Model) 


TX 


Surface  Air  Temperature 
Derived  from  surface  air  and 
potential  temperatures,  boundary 
layer  depth,  upper-level  winds 
extrapolated  to  surface,  air 
density,  drag  coefficient,  gusti- 
ness  factor,  and  empirical  con- 
stants . 


(°C) 


EX 


Surface  Vapor  Pressure 

Derived  from  model's  mixing  ratio, 


SOLARAD   Solar  Radiation 

Calculated  absorption  of  incoming 
short-wave  (solar)  radiation 
(positive  downward) 


(mb) 


(gm-cal/ 

2  u  ^ 
cm  hr ) 
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EHF      Evaporative  Heat  Flux  (gm-cal/ 

Derived  using  air  density,  drag  2  ,  . 

-c-c-  ■   *     4.           -,    j.    j   •  j  cm  hr) 
coefficient,  extrapolated  winds, 

and  mixing  ratios. 

SEHF     Sensible  Plus  Evaporative  Heat  Flux  (gm-cal/ 

SEHF  =  SHF  +  EHF  2  .  , 

cm  hr ) 


SHF      Sensible  Heat  Flux  (gm-cal/ 

Recovered  from  SHF  =  SEHF  -  EHF.  2    . 

Originally  derived  by  FNWC  using        cm     ; 
drag  coefficient,  extrapolated 
winds,  surface  air  temperature,  TX, 
density,  and  constants. 

THF      Total  Heat  Flux  (gm-cal/ 

THF  =  SEHF  -  SOLARAD  +  LW,  where  2  ,  x 

cm   hr ) 
LW  is  the  heating  due  to  long-wave  ' 

(terrestrial)  radiation. 

C.  Marine  Wind  Model  (FNWC) 

WWW     Marine  Wind  Speed  (knots) 

DDWW     Marine  Wind  Direction  (degrees/ 

Both  variables  derived  from  a  10) 

dynamic  balancing  of  surface  wind 
and  sea-level  pressure. 

D.  Spectral  Ocean  Wave  Model  (S .0. W.M . ) (FNWC) 

HW       Significant  Wave  Height  (feet) 

PPW      Primary  Wave  Period  (sec) 

PDW      Primary  Wave  Direction  (degrees/ 

10) 
SPW      Secondary  Wave  Period  (sec) 

SDW      Secondary  Wave  Direction  (degrees/ 

10) 
WCP      Probability  of  White  Caps  (percent) 

E.  Other  Model  Output  Parameters  (FNWC) 

SSTA     Sea  Surface  Temperature  Anomaly         (°C) 
Calculated  anomaly  of  sea-surface 
temperature  from  the  mean  of  the 
day  as  interpolated  from  the 
monthly  mean  values. 
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Derived  Parameters 


u 


CAPU 


Zonal  Wind  Component  (m/sec) 

u  =  -WWW  sin(DDWW  •  10) 

Meridional  Wind  Component  (m/sec) 

v  =  -WWW  cos(DDWW  •  10) 

I  Directional  Wind  Component  (m/sec) 

CAPU  =  -u  •  sin(LNGA)  -  v  •  cos(LNGA), 
(Haltiner,  1971). 


CAPV     J  Directional  Wind  Component 

CAPV  =  u  •  cos(LNGA)  -  v  •  sin(LNGA), 
(Haltiner,  1971),  where  LNGA  = 
-10  -  (I, J  point  longitude). 

THETAX   Potential  Temperature  X 
Derived  using  PS,  TX. 

THETAR   Potential  Temperature  R 
Derived  using  PS,  TAIR. 


(m/sec) 


(°K) 


(°K) 


STABX    Stability  X 

Derived  using  [ THETAX- (THETA  of  T925)] 
/(PS-925).   Value  greater  than  zero 
indicates  absolute  instability. 

STABR    Stability  R 

Derived  using  [THETAR- (THETA  of  T925) ] 
/(PS-925).   Same  value  effect  as  STABX. 

ASTDX    Air-Sea  Temperature  Difference  X 
ASTDX  =  TX  -  TSEA 


(°K/mb) 


(°K/mb) 


(°C) 


ASTDR    Air-Sea  Temperature  Difference  R 
ASTDR  -  TAIR  -  TSEA. 


(°C) 


ADTSEA    Advection  of  TSEA 

Formulae  and  notes  below. 

ADTX     Advection  of  TX 

Formulae  and  notes  below. 

ADTAIR    Advection  of  TAIR. 

Formulae  and  notes  below. 

AASTDX    Advection  of  ASTDX 

Formulae  and  notes  below. 


(°C/ 
Hour) 

(°C/ 
Hour) 

(°C/ 
Hour) 

(°C/ 
Hour) 
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AASTDR   Advection  of  ASTDR 

Formulae  and  notes  below. 

Climatological  Parameters 

CLIMO    National  Climatic  Center  Fog 
Frequency  Climatology 


PERS 


Persistence  of  Fog  Likelihood 


H.   Interactive  Parameters 

CLISEHF  CLIMO  *  SEHF 

CLIV  CLIMO  •  v 

VSEHF  v  •  SEHF 

SEHFSQ  SEHF  •  SEHF 

CLIEHF  CLIMO  *  EHF 

CLIADT  CLIMO  •  ADTAIR 

EHFSQ  EHF  •  EHF 

EHFADT  EHF  •  ADTAIR 


(°C/ 
Hour) 


(percent/ 
100) 

(percent/ 
100) 


(%  gm-cal/ 
100  cm2hr) 

(m  %/ 
100  sec) 

O-00  2m-cal/ 
cm  sec  hr) 

(gm  -cal  / 

2     2 
cm  hr) 

(%  gm-cal/ 

o 
100  cm  hr) 

{%   °C/ 
1Q0  hr) 

(gm  -cal  / 

2     2 
cm  hr) 

(gm-cal   C/ 

cm  hr  ) 


Advection  Formulae  and  Conditions: 

For  the  advection  of  a  quantity  (R)  the  formula,  ADQ  = 
-V'V(Q),  was  used  in  the  following  finite  difference  form: 


RMAP 


AD0- =  -  TsrtCAPn'<*i+i-<Wj  +  capv'^j+i-qj-i)i] 


where        RfiAP  =    ( l+sin(60) )/( l+sin(  latitude) ) 


and 


DM   =    [(2)-(6. 37-10    )«(l+sin(60))]/31. 205 


(31.205  =  grid  mesh  lengths,  pole  to  equator,  on  FNWC's 
I, J  grid). 

In  the  temperature  advection  calculation  for  point  C,  using 
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the  five  grid  points  illustrated  below,  one  or  two  of  the 
points,  namely,  A,  B,  D,  or  E,  may  be  outside  the  study  area. 
In  the  bogusing  method  suggested  by  Mr.  Leo  Clarke,  FNWC, 
when  a  non-center  point  (e.g.,  point  A)  was  judged  to  prob- 
ably produce  a  land  and/or  dissimilar  sea  area  influence  on 
the  resulting  advection,  the  center  point  value  (point  C) 
was  substituted  for  it.   This  "bogusing"  is  necessary  to 
maintain  a  "purely"  marine  characteristic  in  the  resultant 
parameter  value. 

The  set  of  study  area  boundary  grid  points  used  for 
bogusing  (and  some  double  bogusing)  are  depicted  in  Fig.  1. 
The  upper  case  letters  denote  the  points  whose  values  were 
used  for  bogusing.   The  lower  case  letters  mark  the  posi- 
tions of  the  adjacent  points  being  bogused.   The  only  bound- 
ary point  close  enough  to  land  to  give  concern  is  the  one 
near  station  number  32594  (south  tip,  Kamchatka  Peninsula). 
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APPENDIX  C 


Abridged  version  of  internationally  used  weather  code 
figures  and  definitions  for  reporting  present  and  past 
weather  and  low  clouds  in  the  surface  synoptic  report 
(U.  S.  Departments  of  Commerce,  Defense,  and  Transporta- 
tion, 1969). 

Present  Weather 

Code 

Value    Definition 

00-03    Characteristic  change  of  the  state  of  the 
sky  (cloud)  during  the  past  hour. 

04-09    Haze,  dust,  sand,  or  smoke. 

10       Deep  light  fog. 

11-12    Shallow  heavy  fog. 

13-17    Lightning,  thunder,  or  precipitation  within 
sight,  not  reaching  the  ground. 

18-19    Squall(s),  funnel  cloud(s)  during  the  past 
hour. 

20       Drizzle  during  the  past  hour. 

21-23    Rain,  snow,  or  rain  and  snow  during  the  past 
hour. 

24  Freezing  drizzle  during  the  past  hour. 

25-27  Shower(s)  during  the  preceding  hour. 

28  Fog  during  the  past  hour. 

29  Thunderstorm  during  the  past  hour. 

30-39    Duststorm,  sandstorm,  drifting  or  blowing  snow. 

40       Fog  at  distance,  but  not  at  station,  during 
the  past  hour  (visibility  less  than  1  km). 

41-49    Deep  heavy  fog  at  the  time  of  observation 
(visibility  less  than  1  km). 

50-59  Drizzle,  or  drizzle  and  rain. 

60-63  Slight  to  moderate  rain. 

64-65  Heavy  rain. 

66  Slight  freezing  rain. 

67  Moderate  or  heavy  freezing  rain. 

68  Slight  rain  or  drizzle  and  snow. 

69  Moderate  or  heavy  rain  or  drizzle  and  snow. 
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70-79 
80-99 

Past  Weather 


Solid  precipitation  not  in  showers. 

Showery  precipitation  or  precipitation  with 
current  or  recent  thunderstorms. 


0  Cloud  covering  £  or  less  of  sky  throughout  the 
period. 

1  Cloud  covering  more  than  i   of  sky  during  part  of 
the  period. 

2  Cloud  covering  more  than  J  of  sky  throughout  period, 

3  Sandstorm,  or  duststorm,  or  blowing  snow. 

4  Heavy  fog,  thick  haze,  or  smoke. 

5  Drizzle 

6  Rain 

7  Snow,  rain  and  snow. 

8  Shower (s). 

9  Thunderstorm,  with  or  without  precipitation. 

Low  Cloud  Type 

0  No  low  clouds. 

1  Ragged  cumulus  of  fair  weather. 

2  Generally  towering  cumulus. 

3  Cumulonimbus  without  cirriform  or  anvil  tops. 

4  Stratocumulus  formed  by  cumulus  spreading  out. 

5  Stratocumulus  not  formed  by  cumulus  spreading. 

6  Stratus  or  f ractostratus . 

7  Fractostratus  of  bad  weather. 

3    Cumulus  and  stratocumulus,  with  bases  at  different 

levels. 
9    Cumulonimbus  with  cirriform  top. 
/     Low  cloud  obscured. 


Visibility 

90  Less  than  50  m 

91  0-199  m 

92  200-499  m 

93  500  m  -  0.99  km 

94  1  -  1.99  km 

95  2  -  3.99  km 

96  4  -  9.99  km 

97-99  equal  to  or  greater  than  10  km 
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APPENDIX  D 

FOGCAT  I  Categorization  Scheme 
(Van  Orman  and  Renard,  19  77) 

1.   Groupings  and  symbols  used  in  FOGCAT  I  categorization 
scheme . 

Present  Weather  (ww)   Past  Weather  (W)     Low  Cloud  (CL) 


Symbol 

Associated 

Symbol 

Associated 

Symbol 

Associated 

ww  Codes 

W  Codes 

CL  Codes 

41G 

41-49 

4 

4 

6 

6 

10G 

10,28,40 

4,5 

4,5 

5,7 

5,7 

11G 

11,12,20,24 

5G 

0,1,2,5 

B 

/ 

50G 

50-59 

2G 

0,1,2 

* 

1 

any  CL  not 
isted  above 

60G 

60-63,66,68 

* 

any  W  not 
listed  above 

*       any  ww  not 

listed  above 

G  =  Group  B  =  Low  clouds  obscured 

2.   Scheme  for  categorizing  observations  according  to  like- 
lihood of  fog  (FOGCAT  I). 


Major 

Sub- 

Present 

Past 

Low  Cloud 

Category 

Cate 

gory 

Weather 

Weather 

Type 

(ww) 

(W) 

(CL) 

Strong 

al 

41G 

4,5 

6 

Foggers 

a2 

TT 

it 

B 

=  S 

bl 

1! 

2G 

6 

b2 

II 

ii 

B 

b3 

10G 

4,5 

6 

b4 

1! 

ti 

B 

cl 

41G 

* 

6 

c2 

ii 

4,5 

5,7 

c3 

rt 

n 

* 

c4 

10G 

2G 

6 

Foggers 

dl 

41G 

* 

=  F 

d2 

n 

2G 

d3 

50G 

4 

d4 

it 

ii 

d5 

10G 

* 

d6 

ii 

2G 

d7 

it 

4 

B 
5,7 

6 
B 
6 
B 
5,7 
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el 

41G 

* 

e2 

It 

2G 

e3 

50G 

4 

e4 

10G 

* 

e5 

IT 

4 

e6 

11G 

4 

e7 

n 

it 

fl 

41G 

* 

f2 

50G 

4 

f3 

10G 

5G 

f4 

11G 

it 

f5 

!l 

it 

gl 

10G 

* 

g2 

tt 

5G 

g3 

11G 

* 

g4 

it 

4 

g5 

ii 

it 

hi 

ii 

* 

h.2 

IT 

5G 

il 

10G 

* 

i2 

11G 

* 

jl 

it 

5G 

Past /Weak 

kl 

50G 

5G 

Foggers 

k2 

GOG 

4 

=  P 

S-l 

ii 

n 

5,2 

* 

it 

ml 

50G 

* 

m2 

* 

4 

nl 

* 

5G 

ol 

60G 

4 

o2 

* 

4 

Pi 

60G 

4 

p2 

* 

* 

P3 

* 

4 

Maybe 

ql 

50G 

5G 

Foggers 

q2 

ii 

it 

=  M 

q3 

it 

* 

q4 

it 

5G 

q5 

11G 

* 

q6 

60G 

5G 

q7 

it 

* 

rl 

50G 

* 

r2 

it 

* 

r3 

60G 

5G 

r4 

it 

it 

r5 

it 

* 

r6 

it 

* 

5,7 

* 

5,7 
B 

* 

6 

B 

* 

* 

5,7 

6 

B 

5,7 

* 

6 
5,7 

* 

B 
5,7 

* 

5,7 

* 

6 

6 

B 

6 

6 

B 

6 

5,7 

5,7 

* 

6 

* 

B 
5,7 

5,7 

* 

* 

6 
6 
B 

* 

B 

5,7 
B 
5,7 
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Non- 

ul 

60G 

5G 

Foggers 

u2 

it 

* 

=  N 

vl 

* 

5G 

wl 

* 

* 

xl 

* 

* 

yi 

* 

5G 

y2 

* 

5G 

zl 

* 

* 

* 

* 

B 

B 

5,7 

5,7 

* 

Major  fog  categories  and  sub-category  fog  groups  with 
associated  Fog  Probabilities  as  a  function  of 
visibility. 


Visibility  Code  Values 


96- 

-99 

94- 

•95 

90-93 

Fog 

(Good) 

(Fair) 

(Poor) 

Probability 

Major 

Category 

Subcategory 

Groups 

a 

100.0 

a 

b 

96.6 

a 

b 

c 

93.1 

S 

b 

c 

d 

89.7 

c 

d 

e 

86.2 

d 

e 

f 

82.8 

e 

f 

S 

79.3 

f 

S 

h 

75.9 

F 

g 

h 

i 

72.4 

h 

i 

J 

69.0 

i 

J 

65.5 

J 

k 

62.1 

k 

k 

£ 

58.6 

I 

I 

ra 

55.2 

m 

m 

n 

51.7 

P 

n 

n 

o 

43.3 

o 

o 

P 

44.8 

P 

P 

41.4 

q 

q 

q 

37.9 

M 

r 
s 

r 
s 

r 
s 

34.5 

not 

31.0 

used 

t 

t 

t 
u 

27.6 

24.1 

u 

u 

V 

20.7 

V 

V 

w 

17.2 

w 

w 

X 

13.8 

N 

X 

y 

z 

X 

y 

z 

y 

z 

10.3 
6.9 
3.4 
0.0 
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APPENDIX  E 
FOGCAT  II  Categorization  Scheme  (Quinn,  1978) 


Fog 

Present  Weather 

Past  Weather 

Visibility 

Probability 

(ww) 

(W) 

(VIS) 

(%) 

0-9 

_ 

0 

10-12 

— 

- 

100 

13-19 

— 

— 

0 

20 

0-3,6-9 

— 

0 

20 

4,5 

90-95 

35 

20 

4,5 

96-99 

0 

21-24 

0-3,6-9 

- 

15 

21-24 

4,5 

90-95 

60 

21-24 

4,5 

96-99 

15 

25-27 

0-3,6-9 

- 

0 

25-27 

4,5 

90-95 

10 

25-27 

4,5 

96-99 

0 

28 

— 

100 

29-39 

— 

— 

0 

40-49 

— 

— 

100 

50-59 

0-3,6-9 

— 

35 

50-59 

4,5 

90-95 

85 

50-59 

4,5 

96-99 

35 

60-69 

0-3,6-9 

— 

15 

60-69 

4,5 

90-95 

60 

60-69 

4,5 

96-99 

15 

70-79 

— 

0 

80-89 

0-3,6-9 

— 

0 

80-39 

4,5 

90-95 

10 

80-39 

4,5 

96-99 

0 

90-99 

— 

0 

Dash  indicates  that  the  particular  category  was  not  con- 
sidered in  the  assignment  of  a  fog  probability. 
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APPENDIX  F 
Verification  Score  Formulae  (Quinn,  1978) 

EVENT  ESTIMATED 

NO 

C     Total  (T)  =  A+B+C+D 

No.  of  Correct 
D    Forecast  FC  =  A+D 

Range:   ~2BC  „  <  HSS  <  1 

ET  +  CZ  ~ 

(A+B)(A+C)  +  (D+B)(D+C)    „     -       .   , 
where    EX  = — '    m     — -  ,    No.  of  expected 

correct  forecasts  due  to  chance. 
Threat  Score  =  7^  =  J+I+C         Range:   0  <  TS  £  1 

Both  scores  indicate  more  skill  with  larger 
positive  values. 

The  Probability  Score  (PS)  is  from  that  given  by 
Panofsky  and  Brier  (1958)  and  may  be  written  as  (Renard, 
1975): 


YES 

YES 

A 

EVENT 

OBSERVED 

NO 

B 

Heidke 

FC  - 

EX 

Skill  Score 

T  - 

EX 

nl 


P-Score 


-#  i2Z  Pi2  +  E    ^-Pj)2] 

t-i  j=i 


Range:       0    <    PS    <    2 


where 


N  =  Total  number  of  cases 
n   =  Total  number  of  non-events 

p.  =  Associated  probability  value  for  the  non-event 
n,  =  Total  number  of  events 

p .  =  Associated  probability  value  for  the  event 
The  closer  to  zero,  the  greater  the  skill. 
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